Gravitational wave background as a probe of the primordial black hole abundance 
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Formation of primordial black holes (PBHs) requires a large root-mean-square amplitude of den- 
sity fluctuations, which generate second-order tensor perturbations that can be compared with ob- 
servational constraints. We show that pulsar timing data essentially rules out PBHs with 10^~^Mq 
which were previously considered as a candidate of intermediate-mass black hoes and that PBHs 
with mass range IQ^'-'"^^ g may be probed by future space-based laser interferometers. 



Primordial black holes (PBHs) are produced when den- 
sity fluctuations with a large amplitude enters the hori- 
zon in the radiation dominated stage of the early uni- 
verse with their typical mass given by the horizon mass 
at that epoch [l|, Q . PBHs with their mass smaller than 
lO^^g would have been evaporated away by now due to 
the Hawking radiation [3] and a number of cosmologi- 
cal constraints have been imposed on the abundance of 
these light holes by big-bang nucleosynthesis (BBN) Q 
and gamma-ray background [sl] etc. Heavier PBHs can 
play some astrophysical roles today. For example, they 
may serve as the origin of the intermediate-mass black 
holes (IMBHs), which are considered to be the observed 
ultra-luminous X-ray sources, if their mass and abun- 
dance lie in the range Mpbh ^ 10^^« — lO^M© and 
^pbh/^^ - 10-^ - 10-^ respectively [ej, while PBHs 
with mass Mpbh lO^^g - lO^^g (IQ-^^M© - lO-^M©) 
[l, B] and the abundance I^pbh^^ = 0.1 [8] can provide 
astrophysical origin of dark matter (DM) which is yet 
free from the constraint imposed by gravitational lens- 
ing experiments [3, • In order to produce the relevant 
density of PBHs, it is necessary to produce density fluc- 
tuations whose power spectrum has a high peak with an 
amplitude 10~^ — 10~^ on the corresponding scales. A 
number of inflation models has been proposed to realize 
such a spectrum [l2| • 

In these models giving a large amplitude of density 
fluctuations, second-order effects may play important 
roles. For example, they generate non-Gaussianity in 
the statistical distribution of density fluctuation even if 
the linear perturbation is Gaussian. Its effect on PBH 
abundance was recently investigated in single-field infla- 
tion models, but it turned out that non-Gaussian effect is 
negligibly small |13], justifying previous analysis assum- 
ing Gaussianity [ij]. 

Second-order effects also generate tensor fluctuations 
to produce stochastic background of gravitational waves 
(GWs) from scalar-tensor mode coupling 3, 3], whose 
amplitude may well exceed the first-order tensor pertur- 
bation generated by quantum effect during inflation 
in the current set up since the amplitude of density fluc- 
tuations required for a substantial density of PBHs is so 



large. 

In this Letter, we show the GWs induced by scalar fluc- 
tuations as a second-order effect lH, [l6| is a useful probe 



to investigate the abundance of the PBHs. We calculate 
spectrum of these second-order GWs in the case scalar 
fluctuations have a sufficiently large peak to realize for- 
mation of relevant numbers of PBHs. As a natural conse- 
quence we find that the spectrum of GWs has a peak on a 
scale approximately equal to the scale of the peak of the 
scalar fluctuations. We can therefore obtain information 
on the abundance of PBHs with the horizon mass when 
the scale of the peak entered the Hubble radius by observ- 
ing GWs with the frequency corresponding to the same 
comovmes; scale, namely, IQ-^Hz-lQ-^Hz for IMBHs and 
10~^Hz — IHz for dark- matter PBHs. Fortunately, the 
former band is probed by the pulsar timing observations 
0, [l2l] while the latter band can be observed in future 
by the space-based interferometers HqI, [in for dark 
matter PBHs. 

We write the perturbed metric as 

ds^ = a{r]f [-e^^drf + e'^"^ {5ij + hij)dx'dx^] , (1) 

including both scalar perturbations, <l> and ^, and tensor 



perturbation, hij^ which satisfies dih^^ 



with h) 



S'^^hkj' We assume the lowest-order tensor perturbation 
is negligible and consider that generated by the scalar 
mode as a second-order effect. The relevant part of the 
Einstein equation therefore reads 



, d'hi = 2v^p:, (2) 



2nh] 



where a prime denotes differentiation with respect to the 
conformal time, V^j represents the projection operator 
to the transverse, traceless part, and H = a' /a (isl. E^. 
The source term reads 

(3) 

with w = p/p being the equation-of-state parameter of 
the background fluid. In practice, only the radiation 
dominated era is relevant, so we take w = 1/3 hereafter. 



3(1 + ^) 
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We also neglect anisotropic stress, which is expected to 
give only a small correction [16], and set ^ = ^ at lin- 
ear order. In order to calculate the induced GWs up 
to second order, it is sufficient to use the linear scalar 
modes. Hence, we only need to solve the linear evolution 
equation [23], 



V 3 



0, 



(4) 



for the scalar modes where represents a Fourier 
mode. Its non-decaying solution is given by "^kiv) = 
I^/c (77)^^(0) with the transfer function 
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(5) 



For our purpose we assume the following approximate 
form of the power spectrum of the initial fluctuations, 



27r2 



{|*fc(0)|2)=^2^(ln(fc/fcp)), (6) 



where kp and represent the wavenumber of the peak 
and (amplitude)^ x ln(peak width) of the original spec- 
trum, respectively. With this power spectrum the frac- 
tional energy density of the region collapsing into PBHs 
at their formation time is estimated as 



/3(Mpbh) ^ 0.1 exp 



(7) 



where Mpbh is of the order of the horizon mass when 
the comoving scale enters the Hubble radius and 
is the threshold value of PBH formation. Carr [14] takes 



the threshold value of the density contrast to be = 1/3 
corresponding to = 1/2. Analysis based on numerical 
calculation [24| gives a similar but slightly different value 
[l2| . One can express the current value of the density 
parameter of PBHs in terms of /^(Mpbh) as 



^^PBH,o/^' = 2x10^/5(Mpbh) 



V 1036 1 



PBH 



-1/2 



V 10.757 



-1/3 



(8) 



where g^p is the effective number of the relativistic de- 
grees of freedom when hp entered the Hubble radius. 
We define the Fourier modes hk as 



hij{x,T]) 



d^fc 
(27r)3/2 



(9) 

where e^ (/c), e^^- (/c) are polarization tensors which are 

normalized as j ^^j{k)e^-{—k) = 25^^. The Fourier 
transform of the source term (|3]) is also defined similarly. 
We find the source term is constant when /c^t^/a/S ^ 1, 
while it decreases in proportion to 77"^ for kpTj/y/^ ^ 1. 
As a result the production of scalar-induced GWs mostly 
occurs around the time when kp crosses the sound hori- 
zon. Using the Green function method one can easily 
find a formal solution to ([2j) from which we can evaluate 
the density parameter of GWs contributed by a logarith- 
mic interval of the wavenumber around k. It is formally 
expressed as 



(10) 



when /c-mode is well inside the horizon [25|. In the radi- 
ation dominated regime it is explicitly given by 



2 n p 

^Gw{k, ^) = 3 ^Vi / V1V2 sin [k{r] - r]i)] sin [k{r] - 7/2)] Skim^ 



(11) 



where we have defined 



Sk{vuV2)^ / dk / d/i (1 - fi'ffik, \k - klm)f{k. \k - klm)V^{k)V^{\k - k\). (12) 

I 



Here /(/ci, /c2, r/) is a function written in terms of the the spectrum of I^qwI/^^o) today is well described by 
transfer function for the scalar modes as follows. 



f{ki,k2,v) = 2DkMDkM (13) 
- [DkM +n-'D',^irj)][DkAv) +n-'D',^irj)]. 

In the mass range of the PBHs of our interest, creation of 
scalar-induced GWs is terminated well before the equal- 
ity time. After that the energy density of GWs decreases 
in proportion to a~^. As a result the overall feature of 
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for / < /gw = '^fp/V^ where we have used the fre- 
quency / = 27r/c/ao instead of the wavenumber. The 
above expression has the peak value ^^gw(/gw)^^ = 
(16/27)Agw at /gw- 

Note, however, that the actual spectrum of GW cal- 
culated from has a much larger and sharper peak 
at /gw besides the bulk spectrum (fT4j) due to the reso- 
nant amplification between the transfer function and the 
Green function [lH . Such an amplification occurs only if 
the peak width. A, of the primordial scalar fluctuation is 
sufficiently small, A <C kp/2. Since the resonant growth 
of the amplitude depends on the detailed shape of the 
primordial power spectrum around the peak, we do not 
incorporate it to yield a conservative bound on the PBH 
abundance independent of it. 

We now compare our results with observational con- 
straints. For definiteness we identify Mpbh with the 
horizon mass when the peak scale entered the Hub- 
ble radius. This is a reasonable approximation even if 
critical behavior [26] is taken into account Then 
^PBH is related with the peak frequency of GWs as 



logio(MpBH/lg) 
35 30 25 



/gw = 1 X 10-^ Hz 



/MpBH 



-1/2 

Vio.Tsy 



-1/12 



(15) 



The pulsar timing observations are sensitive to GWs with 
/ > 1/T where T is the data span. The 7- year data of 
observation of PSR B 1855+09 gives an upper limit 



^GwU)h^ < 4.8 X 10" 



/ 



4.4 X 10-9Hz 



(16) 



for / > 4.4 X 10-^Hz at 90% confidence level [Til. By 
using this limit, we can constrain the abundance of PBHs 
with mass Mpbh < 1 x 10^^ g = 5 x lO^M©. 

The space-based interferometers are sensitive to GWs 
with 10~^Hz ^ / ^ lOHz, which covers the entire 
mass range of the PBHs which are allowed to be DM, 
lO^^g < MpBH < lO^^g. LISA win have its best sensi- 
tivity Q^G^K^ - 10-^1 at / - lO-^Hz (MpBH ^ lO^^g), 
BBO and the ultimate-DECIGO are planned to have sen- 
sitivity to V^Q^^h? ~ 10~^^ and V^Q^^h? ^ lO"^'', respec- 
tively at / - lO-^Hz (MpBH ^ lO^^g) 

Figure [1] shows the energy density of the induced GWs 
obtained by numerically evaluating (pT]) and tracing its 
subsequent evolution up to the present, whose approxi- 
mate form is given by ([T4j). The left wedge-shaped curve 
represents the case kp = 0.2pc~^ and ^ = 7 x 10~^ cor- 
responding to MpBH = 6 X 1O^M0 and QpBuh'^ = 10"^, 
while the right wedge-shaped curve depicts the case 
kp = 1 X lO^pc"^ and ^ = 9 x 10~^ corresponding to 
MpBH = SxlO^^g and l^PBH^^ = 10"^. Also shown there 
are two envelope curves of the peak value, (16/27)Agw, 
corresponding to I^pbh^^ = 10~^ (red dotted line) and 
^PBH^^ = 10~^ (green broken line) which depend on the 
frequency logarithmically except for the discontinuities 
due to the change of the relativistic degrees of freedom at 
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FIG. 1: Energy density of scalar-induced GWs associated 
with PBH formation together with current pulsar constraint 
(thick solid line segment) and sensitivity of various GW detec- 
tors (convex curves) . Left and right wedge-shaped curves indi- 
cate expected power spectra of GWs from two different peaked 
scalar fluctuations corresponding to (^^pbh/^^, Mpbh, 5'*p) = 
(10-^6 X lO^Mo, 10.75) (left) and (10-\ 3 x lO^^g, 106.75) 
(right), respectively. The red dotted (green broken) line shows 
an envelope curve, ( 16/27) Aqw, corresponding to Qpeh = 
10"^ (10"^) obtained by moving kp and A. 



the QCD phase transition and the electron-positron pair 
annihilation. We have also shown the limit imposed by 
the pulsar timing observation (thick solid line segment) 
and the planned sensitivity of the space-based interfer- 
ometers depicted [28| with the instrumental parameters 
used in [29[ as well as those of LIGO [30]. 

As is seen in the figure the pulsar timing constraint is 
so stringent that one cannot achieve QpBuh^ > 10~^ for 



PBHs with 4 X lO^M© < Mpbh < 5 x lO^M©, ruling 
out the major mass range of IMBHs. On the other hand, 
if pulsar timing experiments should find any nontrivial 
modulation in near future, it might be due to the PBHs 
with mass around lO^M© [l9|. 

It is clear from Fig.[T]that the future space-based inter- 
ferometers can test the feasibility of PBHs as the dom- 
inant constituent of the DM. The ground-based inter- 
ferometers, on the other hand, have good sensitivity at 
/ ~ 10 — lO^Hz [30[. This frequency band corresponds to 
mass scale Mpbh ^ lO^^g — lO^^g. Though the sensitiv- 
ity of LIGO is too low now and in near future to detect 
GWs from the second-order effect associated with PBH 
formation, we could improve the sensitivity by correla- 
tion analysis to reach the desired level to probe PBHs. 
Because the spectrum has a tail to lower frequencies, it 
may be possible to constrain the abundance of the PBHs 
with MpBH < 7 X lO^^g (/gw > 6 x lO^Hz), which have 
evaporated by the present epoch and could contribute to 
cosmic rays. Further study, however, is necessary in order 
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FIG. 2: New constraint on the mass spectrum of PBHs im- 
posed by scalar- generated GWs. Dotted line represents mass 
range to be constrained by future GW detectors. 

to obtain the conclusion because there are astronomical 
sources of GWs in this frequency band, too. 

Figure [2] depicts improved constraint on the PBH frac- 
tion /3(Mpbh) where dotted region denotes the mass 
range to be constrained by future laser interferometers. 

In summary, we have calculated the spectrum of the 
stochastic gravitational wave background generated as 
a second-order effect from scalar perturbations which 
have a spectrum with a high peak to realize formation 
of appreciable numbers of PBHs. As a result we have 
found that PBHs with their mass corresponding to that 
of IMBHs are already being ruled out because the ampli- 
tude of the associated GWs exceeds the limit imposed by 
the pulsar timing. We have also found that if PBHs with 
mass 10^^ — lO^^g are dominant constituents of DM, we 
can easily detect the relevant GWs by future space-based 
laser interferometers. Thus gravitational waves are new 
and powerful probe the mass spectrum of PBHs. 

We thank K. Ichiki for useful comments. This work 
was supported in part by JSPS Grant-in- Aid for Scientific 
Research No. 19340054(JY) and by Global COE Program 
"the Physical Sciences Frontier" , MEXT, Japan. 



[1] S. Hawking, Mon. Not. Roy. Astron. Soc. 152, 75 (1971); 
B. J. Carr and S. W. Hawking, Mon. Not. Roy. Astron. 
Soc. 168, 399 (1974). 

[2] B. J. Carr, In Proc. 59th Yamada Conference "Inflat- 
ing Horizons of Particle Astrophysics and Cosmology", 
(Universal Academy Press, Tokyo, 2005) p 129. 

[3] S. W. Hawking, Nature 248, 30 (1974). 

[4] Ya. B. Zeldovich, A. A. Starobinsky, M. Yu. Khlopov 
and V. M. Chechetkin. Pis'ma Astron. Zh. 3, 208 (1977) 
[English translation: Sov. Astron. Lett. 3, 110 (1977)]; 



K. Kohri and J. Yokoyama, Phys. Rev. D 61, 023501 
(2000). 

[5] D. N. Page and S. W. Hawking, Astrophys. J. 206, 1 

(1976); J. H. MacGibbon and B. J. Carr, Astrophys. J. 

371, 447 (1991). 
[6] T. Kawaguchi, M. Kawasaki, T. Takayama, M. Yam- 

aguchi and J. Yokoyama, Mon. Not. Roy. Astron. Soc. 

388, 1426 (2008). 
[7] M. A. Abramowicz, J. K. Becker, P. L. Biermann, 

A. Garzilli, F. Johansson and L. Qian, larXiv:0810.3140l 

[astro-ph]. 

[8] E. Komatsu et al. arXiv:0803.0547 [astro-ph]. 
[9] C. Alcock et al [MACHO Collaboration and EROS Col- 
laboration], Astrophys. J. Lett. 499, L9 (1998); P. Tis- 
serand et al. [EROS-2 Collaboration], Astron. Astrophys. 
469, 387 (2007). 

[10] G. F. Marani, R. J. Nemiroff, J. P. Norris, K. Hurley and 
J. T. Bonnell, Astrophys. J. Lett. 512, L13 (1999). 

[11] P. Ivanov, P. Naselsky and I. Novikov, Phys. Rev. 
D 50, 7173 (1994); J. Garcia-Bellido, A. D. Linde 
and D. Wands, Phys. Rev. D 54, 6040 (1996); 
J. Yokoyama, Astron. Astrophys. 318, 673 (1997); 
M. Kawasaki and T. Yanagida, Phys. Rev. D 59, 043512 
(1999); M. Kawasaki, T. Takayama, M. Yamaguchi and 
J. Yokoyama, Phys. Rev. D 74, 043525 (2006). 

[12] J. Yokoyama, Phys. Rev. D 58, 083510 (1998); Phys. 
Rept. 307, 133 (1998); Prog. Theor. Phys. Suppl. 136, 
338 (1999). 

[13] R. Saito, J. Yokoyama and R. Nagata , JCAP 0806, 024 

(2008); J. C. Hidalgo, arXiv: 0708:3875 [astro-ph]. 
[14] B. J. Carr, Astrophys. J. 201, 1 (1975). 
[15] K. N. Ananda, C. Clarkson and D. Wands, Phys. Rev. D 

75, 123518 (2007). 
[16] D. Baumann, P. J. Steinhardt, K. Takahashi and 

K. Ichiki, Phys. Rev. D 76, 084019 (2007). 
[17] A.A. Starobinsky, JETP Lett. 30, 682 (1979); V.A. 

Rubakov, M.V. Sazhin, and A.V. Veryaskin, Phys. Lett. 

115B, 189(1982); L.F. Abbott and M.B. Wise, Nucl. 

Phys. B 244, 541 (1984). 
[18] S. E. Thorsett and R. J. Dewey, Phys. Rev. D 53, 3468 

(1996). 

[19] F. A. J enet et al, Astrophys. J. 653, 1571 (2006). 
[20] Ihttp777lisa.nasa.gov/ 

[21] N. Seto, S. Kawamura and T. Nakamura, Phys. Rev. 
Lett. 87, 221103 (2001); S. Kawamura et al, J. Phys. 
Conf. Ser. 120, 032004 (2008). 

[22] S. Phinney et al., The Big Bang Observer: Direct De- 
tection of Gravitational Waves from the Birth of the 
Universe to the Present, NASA Mission Concept Study, 
2004. 

[23] V. F. Mukhanov, H. A. Feldman and R. H. Branden- 
berger, Phys. Rept. 215, 203 (1992); H. Kodama and 
M. Sasaki, Prog. Theor. Phys. Suppl. 78, 1 (1984). 

[24] M. Shibata and M. Sasaki, Phys. Rev. D 60, 084002 
(1999); I. Musco, J. C. Miller and L. Rezzolla, Class. 
Quant. Grav. 22, 1405 (2005). 

[25] M. Maggiore, Phys. Rept. 331, 283 (2000). 

[26] J. C. Niemeyer and K. Jedamzik, Phys. Rev. Lett. 80, 
5481 (1998). 

[27] J. Yokoyama, Phys. Rev. D 58, 107502 (1998). 

[28] http://www.srl. caltech.edu/^ shane/sensitivity/ 

[29] H. Kudoh, A. Taruya, T. Hiramatsu and Y. Himemoto, 

Phys. Rev. D 73, 064006 (2006). 
[30] jhttp://www.ligo.caltech.edu/^ 



